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ABSTRACT 


This  report  concerns  a  transform  based  on  a  set  of  two- 
dimensional  Haar-like  functions.  Series  expansions  of  two-dimensional 
functions  in  terms  of  this  set  have  convergence  properties  that  are 
analogous  to  those  of  Haar  series.  Following  a  brief  review  of  Haar 
functions,  the  two-dimensional  functions  are  defined  and  their  prop¬ 
erties  discussed.  Analog  and  digital  transforms  based  on  these  func¬ 
tions  are  then  considered;  Possible  applications  are  image  coding 
and  pattern  recognition. 
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A  TWO-DIMENSIONAL  HAAR-UKE  TRANSFORM 


1.  INTRODUCTION 

In  two-dimensional  image  processing,  use  of  the  two-dimensional  Fourier  transform  is 
being  supplemented  by  the  Walsh-Haramard  transform  (1).  Haar  functions  (2-4)  are  closely 
related  to  Walsh  functions  (5)  and  have  properties  that  suggest  they  may  be  useful  in  ap¬ 
plications  such  as  multiplexing,  data  transmission,  and  pattern  recognition  (6).  Recently, 
a  set  of  two-dimensional  functions  with  properties  analogous  to  Haar  functions  has  been 
introduced  (7).  These  two-dimensional  functions  may  be  well  suited  to  image  coding  and 
related  applications. 

Haar  Functions  are  briefly  reviewed  in  Section  2.  In  Section  3  we  describe  the  two- 
dimensional  functions.  Analog  and  digital  transforms  are  discussed  in  Section  4.  Compu¬ 
tational  and  memory  requirements  are  analyzed.  Applications  are  discussed  in  Section  5. 
Before  continuing,  we  introduce  some  definitions  that  will  be  used  throughout,  this  report. 

Definition.  Let  the  interval  ((n  -  l)/2n_1,  n/2n-1 )  be  denoted  by  i”,  where  m  =  1, 
2, ....  2n_1.  Thus  the  interval  [0, 1]  is  subdivided  into  2a_1  equal  parts  labeled  i„,  i2,  ..., 
In  • 


Definition.  Let  the  sauare  defined  by  (£  -  l)/2n_1  <  x  <  £/2n_1  and  (k  -  < 

y  <  k/2n_1  be  denoted  by  a*8,  where  k,  8  =1,  2, ...,  2n-1.  The  square  is  the  inter¬ 
section  of  i8  along  the  x  axis  and  of  i„  along  the  y  axis.  The  unit  square  will  be  denoted 
!•  When  the  boundaries  of  the  square  are  to  be  included,  one  writes  [a„C]  instead  of 

<vu  . 


Definition.  If  x  E  [0, 1]  for  some  k  and  M  satisfies  x  *  k/2M,  where  k  *  0, 1, 

2, ...,  2m,  then  is  called  a  binary  rational.  If  x  does  not  satisfy  x  =  k/2M  for  any  M, 
then  it  is  called  a  binary  irrational.  A  point  (x,  y)  E  is  called  a  binary  rational  if  either 
x  or  y  is  a  binary  rational.  It  is  called  a  binary  irrational  if  both  x  and  y  are  binary  irra¬ 
tionals.  A  binary  rational  line  segment  is  c  set  of  colinear,  binary-rational  points  with  seg¬ 
ment  endpoints  that  are  binary  rational  in  both  x  and  y. 

Definition.  Let  f(x,  y)  e  L2^];  then 


f(o J8)  =  22”"2 


n 


f(x,y)dxdy 


is  the  mean  value  of  f  over  the  square  . 


(1) 
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\  DEFlNmoN.  If  t,  g  €  L2  [0, 1] ,  then  let 

\  (n>) s 

i 
i 

t 
! 

|  If  f,  g  €E  L2[<hJ,  fnen  tef 

|  (/I*)  =  j*  f(x,  y  )g(x,  y)  dxdy.  (3) 


2.  BRIEF  REVIEW  OF  HAAR  FUNCTIONS 

The  Haar  orthonormal  sequence  is  defined  on  the  closed  interval  [0, 1]  and  is  com¬ 
posed  of  functions  labeled  by  two  indices: 

MT>  ■  . (4) 

The  functions  are  defined  as  follows: 


f  f(x)g(x)dx.  (2) 
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At  points  of  discontinuity,  the  Haar  functions  are  defined  to  be  the  average  of  the  limits 
approached  on  the  two  sides  of  the  discontinuity.  The  first  few  Haar  functions  are  shown 
in  Fig.  1.  We  note  that  <^{x)  is  nonzero  over  the  interval  i™.  The  Haar  functions  are  a 
complete  orthonormal  basis  of  L2[ 0, 1]. 


Fig.  1— Tht  first  tight  Haar  functions 


Any  function  can  be  expressed  as  an  infinite  series  in  terms  of  Haar  functions: 

00  Jrt-l 

/(*)  "  C0  +  Zj  C"  ' *6) 

n»l  m»l 

where  c™  *  (fl^n  ).  For  the  purposes  of  this  report,  convergence  is  best  discussed  by 
means  of  the  partial  sums 
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,  N  2n~1 

i  sN(X) = co +  2]  J2  c”^«(x)’  (6) 

n-  1  m»l 

t 

1  which  contain  2N  terms.  Also  SN+1  contains  2N  more  terms  than  SN,  namely  all  Haar 

I  functions  with  the  subscript  N  +  1.  More  general  partial  sums  are  discussed  in  (2)  and 

j  (4). 

r 

|  For  continuous  functions,  the  sequence  of  partial  sums  SN  is  uniformly  convergent  to 

>  the  given  function.  For  discontinuous  functions,  {&#}  will  still  converge  uniformly  to 

(1/2 )[f<x+)  +  /(*-)] ,  provided  all  discontinuities  are  at  binary-rational  points.  This  con¬ 
vergence  property  for  discontinuous  functions  derives  from  the  fact  that  all  Haar-function 
discontinuities  are  at  binary -rational  points.  For  functions  with  discontinuities  at  binary- 
irrational  points,  SN(x),  though  no  longer  uniformly  convergent,  is  still  pointwise  conver¬ 
gent  everywhere  except  at  the  binary-irrational  discontinuities. 

Several  aspects  of  the  potential  utility  of  Haar  functions  derive  from  an  important 
property  of  the  partial  sum.  In  the  expansion  of  f(x),  the  Nth  Haar  partial  sum  SN(x) 
is  a  step  function  with  2N  equal-length  steps.  The  value  of  S^(x)  on  each  step  is  the 
mean  value  of  f(x)  in  the  interval  covered  by  the  step.  The  value  of  SN{x)  at  a  discon¬ 
tinuity  between  adjacent  steps  is  halfway  between  the  adjacent  steps.  Stated  differently, 
SN(x)  is  the  step  function  of  2N  steps  which  is  the  best  approximation  to  f(x)  in  the  mean- 
square -error  sense.  The  effect  of  additional  terms  is  simple  and  intuitive,  unlike  the  effect 
of  additional  terms  when  the  function  is  expanded  as  a  trigonometric  Fourier  series  (or  as 
a  aeries  in  terms  of  other  continuous  bases  of  L2[ 0, 1]).  We  note  that  it  follows  from 
this  mean-value  property  that  if  f(x)  is  constant  in  the  interval  covered  by  any  step,  then 
Sff(x)  ■  f(x)  exactly  on  this  step. 

The  coefficients  in  the  Haar  Series  also  have  a  simple  relationship  with  the  mean 
value  of  f(x)  over  subintervals  of  [0, 1] .  The  coefficient  Cg1  is  proportional  to  the  differ¬ 
ence  of  two  adjacent  steps  of  Sg(x),  namely  the  steps  on  either  side  of  x  *  (2m  - 1)/2®. 
This  property  leads  to  the  modest  computational  requirements  of  the  Haar  transform' (6). 


3.  TWO-DIMENSIONAL  HAAR-LIKE  FUNCTIONS 

The  attractive  qualities  of  Haar  functions  derive  from  two  related  properties.  First, 
in  the  expansion  of  f(x),  SN(x)  is  a  step  function  of  2N  equal-length  steps.  The  value  of 
SN  on  each  step  is  the  mean  value  of  /  in  the  interval  covered  by  the  step.  Second,  the 
expansion  coefficients  are  proportional  to  the  difference  in  the  mean  value  of  f(x)  over 
adjacent  subintervals  of  [0, 1].  As  discussed  in  (6),  these  properties  may  be  useful  in  ap¬ 
plications  such  as  data  coding  and  pattern  recognition.  Whatever  the  application,  Haar 
functions  require  a  one-dimensional  data  stream.  If  the  data  is  intrinsically  two  dimen¬ 
sional,  it  seems  more  appropriate  to  use  a  two-dimensional  set  of  functions  with  analogous 
properties. 

We  therefore  look  for  a  set  or  orthonormal  functions  that  have  the  following  proper¬ 
ties.  When  a  function  on  the  unit  square  f(x,  y )  is  expanded  in  termr  of  these  functions, 

\ 
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the  Nth  partial  sum  Pw(x,  y)  is  a  step  function  of  2W  square  steps  each  of  area  1/22JV. 
The  value  of  P^(x,  y)  on  any  step  is  the  mean  value  of  f(x,  y)  over  the  area  covered  by 
the  step.  The  expansion  coefficients  are  proportional  to  the  difference  in  the  mean  value 
of  f(x,  y)  over  adjacent  subareas  of  the  unit  square. 

In  analogy  with  the  double  trigonometric  Fourier  Series,  we  might  expect  that  such 
a  set  is  provided  by  all  products  of  the  form 

y)  ■  *5?(*)*J<y).  (7) 

The  set  }  is  the  basis  of  the  two-dimensional  Haar  transform  used  in  Ref.  8.  In  fact, 
no  subset  of  these  products  defines  an  adequate  two-dimensional  sequence,  in  the  sense 
described  above.  Instead,  we  find  that  the  desired  sequence  is  obtained  by  defining  three 
sets  of  functions:  S^(x,  y),  VlJ(x,  y),  and  J ffff(x,  y).  The  first  is  a  subset  of  the  products 
defined  in  Eq.  (7): 

SQ(x,  y)  «  <fi0(y)<fio(x)  (8) 

siVy)*^4(y)#)  O) 

where  (x,  y)  €  and  n  >  0.  The  first  few  are  shown  in  Fig.  2.  The  symbol  S  was 
chosen  because  the  functions  have  a  saddle  shape. 

Introducing  the  block  functions 

2<n-iy2,  x<=tm 

<(*)  -  (10) 

0,  elsewhere 

we  define 

H%(X,  y)  -  ¥>n(y )<*»(*)  (11) 

vU(x,y)  -oi(y)ri(*)  (12) 

where  (x,  y)  €  a*/  and  n  >  0.  The  first  few  are  shown  in  Fig.  3.  The  symbol  H  was 
chosen  because  the  rectangular  nonzero  areas  are  horizontally  oriented.  The  first  few 
Vf/  are  shown  in  Fig.  4.  The  symbol  V  was  chosen  because  the  rectangular  areas  are 
vertically  oriented.  The  form  of  the  definitions  was  chosen  so  that  the  superscripts  locate 
the  nonzero  square  &{/  according  to  the  familiar  row-column  matrix  notation. 

On  the  four  borders  of  of/,  all  sj/,  H%,  and  V*J  are  defined  as  the  average  of  the 
limits  approached  on  the  two  rides  of  the  border.  For  example,  if  has  no  borders  in 
common  with  the  unit  square  Oj,  then  the  functions  take  one  of  the  values  0,  ±2R~2.  On 
a  common  border  of  and  they  take  one  of  the  values  0,  ±2B-1. 


Fig.  2-The  first  few  saddle  function*  S$(x,  y).  The 
shaded  are».  H*ve  the  value  zero.  Other  areas  have  value* 
of  plus  or  minus  the  constant  written  outside  the  square. 


.8  S 


■tmmi 
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Fig.  4— The  first  few  vertical  function*  Vtfx,  y) 

As  shown  in  Ref.  7,  the  set  {<$«(*,  y),  Vfj(xt  y ) }  is  orthonormal  and  complete 
I/2[<vi] .  Series  expansions  in  terms  of  this  set  have  convergence  properties  analogous 
those  of  one-dimensional  Haar  series.  The  Nth  partial  sum  is  given  by 

n  a"-1  a""1  • 

Pn(*>  y)  *  Vo(*.y)  +  Z]  £  L  [°",s"i(x'  *y*  c"  v«  <*•  »>]  ■ 


n«l  i-1  /-I 


where 


fl0  "  (f\S0) 

4  -  (n*> 

»?--  (/I*?) 

'!!  -  (/iv?) 

Ptf(x,  y) is  a  step  function  of  2™  square  steps  each  of  area  1/2^.  The  value  of 
Pfi(xt  y)  on  any  step  is  the  mean  value  of  f{ x,  y)  over  the  area  covered  by  the  step.  If 
f(x,  y)  is  continuous  or  has  a  finite  number  of  discontinuities  along  binary-rational  line 
segments,  then  PN(x,  y)  converges  uniformly.  If  f(x,  y)  has  a  finite  number  of  discontinui¬ 
ties  along  binary-irrational  line  segments,  then  P^ix,  y)  converges  pointwise,  except  along  the 


•’-*»»  V'S'%  V  *  „  i. 


tv-**** 
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discontinuities.  The  accuracy  of  the  estimate  Pw(x,  y)  in  the  neighborhood  of  some  point  is 
proportional  to  the  gradient  at  that  point.  Also  the  accuracy  of  the  estimate  is  bounded 
by  a  value  which  is  half  that  which  bounds  the  estimate  Ps.  The  functions  S^,  HlJ,  and  Vff 
are  nonzero  only  over  the  square  Each  of  the  coefficients  aj{,  bff,  and  cff  can  be  written 
in  terms  of  the  mean  value  of  f(x,  y )  over  the  four  quadrants  of  (Fig.  5); 

«r  ■  2"‘‘)  -  «** s”)  -  ft*2”-1)  +  ft*2™)]  <i4» 

»r  ■  £  [ft*1  ♦  ft*1  2”)  -  ft*2”’1)  -  ft*2”)]  (14b) 

**  *  i  »  v  ‘  m  / 

„ nm  m  ^  [t/ .  2n-l  2m~l\  -  2n-l  2m  v.  7/.’Sta.2»n-l .  7,.2n  2m  J 

ci  )  ~T(<^+i  )  +  r(«S+1  )  “  )J  (14c) 


m-l  2m-l  m 

•y^,-  a**1  '.a!  24-1 


■4 

«f  '  ' 

,2rrl  2m-i 
*!♦» 

• '  ts-l  2m 
'<♦1 

,2b  2*-| 

''I'M 

.2  tsi 

Fig.  6— The  aqutre  6gm  and  its  four  quarters 

4.  TRANSFORMS 
4.1  Analog 

A  straightforward  means  of  obtaining  the  coefficients  0jm ,  b% m ,  and  Cjm  is  suggested 
by  Eqs.  (14a)-(14c),  as  well  as  by  Figs.  2  -  4.  In  the  transform  of  some  function  f(x,  y),  each 
coefficient  is  expressed  as  sums  and  differences  of  the  mean  value  of  f  over  different  areas. 

If  we  let  /(<$)  represent  the  integral  of  f(x,  y)  over  the  square  then  Eq.  (14a)  may  be 
rewritten  as 

«r  -  [k*1  2m-1)  -  K*'  2">  -  K*2"-1)  +  ft*2™)],  aw 

Equations  (14b)  and  (14c)  may  be  rewritten  similarly. 
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Dropping  for  the  moment  the  factor  26"1,  each  coefficient  in  the  transform  of  a 
photograph  can  be  obtained  as  follows.  A  light  beam  is  passed  through  the  negative  and 
then  split.  Each  of  the  two  resulting  beams  is  interrupted  by  a  mask.  Light  not  blocked 
by  the  maak  is  collected  at  a  suitable  integrating  detector,  perhaps  a  photomultiplier.  The 
two  detector  outputs  are  subtracted,  and  the  result,  one  coefficient  of  the  transform,  is 
converted  to  digital  form.  Il'  required,  the  factor  2s”1  may  be  inserted  by  means  of  an 
(S  -  l)-bit  left  shift.  As  an  example,  masks  for  obtaining  a^1,  b^1,  and  c^1  are  shown 
in  Fig.  6. 


An  alternative  method  uses  only  one  beam  and  detector:  the  negative  is  sequentially 
exposed  to  the  two  masks,  the  detector  outputs  are  stored  digitally,  and  the  subtraction 
is  performed  digitally.  Although  slower,  this  method  would  be  cheaper  and  would  avoid 
balance  problems. 


4.2  Digital 

We  now  cousider  the  case  where  the  photograph  has  been  encoded  as  «  *  2^  *  4N 
samples  and  the  transform  is  to  be  obtained  digitally.  This  data  might  be  obtained  by 
digitizing  the  output  of  a  flying-spot  scanner.  Bearing  in  mind  that  each  sample  is  usually 
the  average  value  over  the  smallest  resolution  cell,  and  in  analogy  with  the  modified  Haar 
transform  discussed  in  Ref.  6,  we  define  a  modified  transform  as  follows: 
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$nm  =  2W4+lanm 

(16a) 

ym  =  22N-«+l6»m 

(16b) 

gmn  =  22N-t+lcnm 

(16c) 

When  £  =  N,  Eq.  (16a)  corresponds  to  dropping  the  leading  factor  in  Eq.  (14a).  When 
£  <  N,  it  corresponds  to  replacing  the  leading  factor  in  Eq.  (15)  with  the  constant  2^, 
which  is  needed  to  compensate  for  the  fact  that  the  initial  samples  are  averages  over  the 

smallest  area  of  resolution  1/2^.  ... 

*  >*-'  ’ 

For  illustration  we  discuss  the  case  N  ‘  We  takevthe  n  =  42  =  16-point  transform 

of  data  given  by  a  four-by-fdur  miatrix  x^.  The  brute-force  calculation  then  proceeds  as 
follows:  .  . • 

a0  -  xn  +  x12  +  x18  +  x14  +  *21  +  * 22  +  *23  +  *24 

+  +  *32  +  *33  +  *34  +  *41  +  *42  +  *43  +  *44 

fil1  =  *11  +  *12*  -  *13  “  *14  +  *21'  +  *22  ~  *23  -  *24 

-  *3i  ”  *82  +  *33  +  *34  ~  *41  "  *42  +  *43  +  *44 

&J1  -  XU  +  *12  +  *13  4  *14  +  *21  +  *22  +  *23  +  *24 

"  *31  ”  *32  "  *83  "  *84  *  *41  ”  *42  “  *43  “  *44 

etc.  This  calculation  requires  112 'additions.  In  general,  Sn  log4  n  (1  +  3  log4)  additions 
are  required  for  an  n  ■  4^-point  trarisifbrm. 

As  with  the  Fourier,  Walsh,  andHaar  transforms,  a  fast  transform  results  from  the 
proper  grouping  of  terms.  As  in  the  fast  modified  Haar  transform,  sums  and  differences 
are  calculated  at  each  stage: 


“1 

-  *11 

+  *12 

Vt  • 

1  *11 

-  *12 

«3 

•  *21 

+  *22 

vs  " 

*21 

-  *22 

“6 

“  *31 

+  *32 

y5  " 

1  *31 

"  *32 

“7 

*  *51 

+  *42 

V7  - 

1  *41 

-  *42 

“  V1 

-  vs 

hf 

-  «i 

-  «3 

*f 

“  v2 

-  v4 

K2 

-  «2 

“  “4 

a? 

“  Vg 

-  Vij 

bf 

"  «6 

-  u7 

if 

-  Vg- 

■  w8 

bf 

"  “6 

-  Ug 

« 2  “  *13  +  *14 
“4  “  *23  +  *24 
«6  -  *33  +  *34 
“8  *  *43  +  *44 

5“  l>i  +  V3 
if  -  V2  +  VA 
af1  -  i»6  +  v7 

if  «  i/6  +  Vg 


V2  -  *13  -  *14 
t>4  -  *23  “  *24 
H  “  *33  -  *34 
<>8  “  *43  -  *44 

«>1  ■  «1  +  “3 
W2  M  U2  +  U4 
Wg  •  U5  +  u7 
W4  -  U6  +  Ug 


wwW  t'lwwwuww —————— — - - -r- ^ -~ - - - — - - - - ... - - — — **** — +  w****-.^'  %'**y?*^ **' 


yi 

y2 


*11 

°1 


Wi  +  U>2 
u>3  +  it;  4 
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*1  T  “>1  -  u,2 
z2  =  u>3  -  w4 


Zi  -  z2  &}1  =  yx  - 


y2 


gj1  -  Zi  +  z2 
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50  ■  yi  +  y2 


This  fast  transform  requires  only  40  additions.  The  general  requirement  is  (8/3)(n  -  1). 


Since  it  is  similar  to  the  modified  Haar  transform,  the  modified  transform  just  presented 
is  appropriate  for  applications  such  as  pattern  recognition.  No  information  is  lost  by 
modifying  the  leading  factors  in  Eqs.  (14a)-(14c)  and  (15),  since  the  correct  factor,  which 
is  given  by  the  identity  of  the  coefficient,  can  always  be  reinserted.  We  should  be  careful, 
however,  in  using  the  mo-  ified  transform  in  applications  that  involve  transmitting  it  over 
a  noisy  channel.  Depending  on  the  coding  technique  and  on  the  nature  of  the  channel, 
use  of  the  modified  transform  can  result  in  unequal  expected  error  rates  for  different  co¬ 
efficients.  This  is  because  unequal  energies  may  be  used  to  transmit  different  coefficients. 

The  complete  two-dimensional  transform  is  obtained  by  multiplying  the  results  of  the 
modified  transform  by  the  correct  factor  2^"2Ar“1  (see  Eqs.  (16a)-(16c).  Assuming  a 
binary  representation,  this  can  be  done  by  a  (2N  -  2  +  l)-bit  right  shift.  One  multibit 
shift  is  needed  for  each  coefficient,  so  that  n  multibit  shifts  are  needed  in  the  n  =  4W- 
point  transform.  The  computational  requirements  for  the  modified,  modified  fast,  and 
complete  fast  two-dimensional  transforms  are  summarized  in  Table  1. 


Table  1 

Computational  Requirements  for  n  =  4^ -point, 
Two-Dimensional  Transforms 


Transform 

Additions 

Multibit 

Shifts 

Modified 

3n  log4  n 

• 

Modified  Fast 

!<»-!> 

- 

Complete  Fast 

!<"-» 

n 

In  the  fast  Fourier  and  Walsh  transforms,  the  average  number  of  operations  per  point 
increases  as  log2  n,  so  that  the  arithmetic  speed  requirement  is  a  function  of  both  the  data 
rate  and  the  transform  size.  In  the  fast  Haar  transform,  the  average  number  of  operations 
par  point  is  independent  of  the  transform  size.  The  same  is  true  of  the  two-dimensional 
Haar-like  transform  just  introduced.  The  arithmetic  speed  requirement  is  a  function  of  the 
data  rate  alone.  The  only  limitation  on  transform  size  is  that  imposed  by  the  amount  of 
available  storage. 


“•*  *  «» 
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If  the  n  *  4n  samples  are  located  in  memory  prior  to  the  transform,  this  is  sufficient 
space  in  which  to  complete  the  transform.  If  the  samples  are  accepted  one  at  a  time  from 
an  external  source  end  if  transformed  coefficients  can  be  put  out  immediately  after  cal¬ 
culation.  then  the  memory  requirement  is  reduced.  The  extent  of  this  reduction  depend).' 
on  the  order  in  which  the  samples  are  made  available  for  coir— itation.  If  the  samples 
derive  from  a  raster  scan  pattern  in  which  a  full  row  of  data  i*  supplied  before  flyback, 
then  2\/n-  2  memory  locations  are  required  for  the  transform.  If  an  optimum  scan  pat¬ 
tern  is  used,  the  memory  requirement  drops  to  3  log4  n  locations.  Optimum  scan  pat¬ 
terns  encode  successively  larger  squares,  starting  at  one  comer  of  the  image.  For  example, 
the  following  is  an  optimum  scan  pattern  for  the  16-point  transform  discussed  previously: 

*11»  *12*  *21*  *22*  *13*  *14*  *23*  *24*  *31*  *32*  *41*  *42*  *33*  *34*  *43*  *44* 

The  calculation  then  takes  place  in  the  following  order: 

**1.  *1.  «3.  &21.  «>1.  VS,  Sjj1,  cl1, 

«2*  v2>  “4*  &22*  w2»  *>4*  S22*  *22* 

U5,  t>5*  “7*  &21*  w3>  y7*  4*.  62l* 

«6*  y6*  “8*  &21*  »4»  y8»  S22*  *22* 

yi»  y2*  M1*  ao»  *i*  *2.  ^l1*  *?• 

Here  six  storage  locations  are  needed,  sinoe  at  one  stage  of  the  calculation  wlt  w2,  w3, 
u6,  v(,  and  u8  must  be  retained. 


5.  APPLICATIONS 


The  application  of  the  two-dimensional  Haar-like  functions  to  two-dimensional  data 
is  analogous  to  the  application  of  Haar  functions  to  one-dimensional  data. 


One  example  is  image  coding.  Information  contained  in  an  image  is  often  transmitted 
by  means  of  the  coefficients  of  an  expansion  in  terms  of  some  set  of  basis  functions.  If 
convergence  is  rapid  or  if  many  coefficients  are  zero,  we  can  often  reduce  the  transmission 
bandwidth  from  that  required  to  send  the  image  itself  (1). 


In  the  two-dimensional  Haar-like  transform,  a  given  coefficient  a%,  6^ ,  or  c^  contains 

The  full  set  of  4N  coefficients  may  be  said  to  con- 

li  ** 


information  only  from  the  square 
tain  a  mixture  of  local  and  global  information.  All  points  contribute  to  a0,  a}1,  and 
cj1;  one  quarter  of  the  points  contribute  to  aj1;  etc.  In  general,  each  point  contributes 
to  3N  of  the  4N  coefficients.  One  result  of  this  is  that  there  is  less  immunity  to  channel 
errors  than  in  the  case  of  a  transform  in  which  all  points  contribute  to  every  coefficient, 
such  as  the  Fourier  transform.  On  the  other  hand,  greater  bandwidth  reduction  may  be 
possible  if  only  the  significantly  nonzero  coefficients  are  transmitted. ,  For  example,  assum¬ 
ing  a  4n  point  transform,  if  a  function  is  constant  over  the  square  ,  where  Af  <  N,  then 
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coefficients  are  identically  zero. 

The  potential  for  bandwidth  reduction  is  summarized  by  the  statement  that  coeffi¬ 
cients  of  the  transform  are  proportional  to  the  differences  in  the  mean  value  of  fix,  y) 
over  acfjacent  squares  of  area  1/4** ,  where  k  -  0, 1,  2,  ...,N  (see  Eqs.  (14a)-(14c)).  The 
two-dimensional  Haar-like  transform  should  be  particularly  appropriate  for  the  transmission 
of  images  that  have  relatively  large  areas  of  constant  or  slowly  changing  tone. 

Equations  (14a)-(14c)  also  suggest  that  the  transform  may  be  useful  for  edge  detec¬ 
tion  in  two-dimensional  images.  The  coefficients  are  sensitive  to  horizontal  edges;  the 
coefficients  are  sensitive  to  vertical  edges;  and  the  coefficients  are  sensitive  to  saddle 
points.  Stated  differently,  {ffjjj}  is  a  set  of  horizontal  edge  detectors;  {V^}  is  a  set  of 
vertical  edge  detectors;  and  {S‘i)  is  a  set  of  saddle  detectors.  Because  edge  detection  by 
differencing  is  sensitive  to  noise,  if  the  (sj£,  ff|£,  V^j}  set  is  useful,  it  is  likely  to  be  so 
only  as  one  part  of  an  edge-detection  process  (9, 10). 
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